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We employ the extended 1 + 3 orthonormal frame formalism for fluid spacetime geometries 
(Ai, g, u), which contains the Bianchi field equations for the Weyl curvature, to derive a 44-D 
evolution system of first-order symmetric hyperbolic form for a set of geometrically defined dynamical 
field variables. Describing the matter source fields phenomenologically in terms of a barotropic 
perfect fluid, the propagation velocities v (with respect to matter-comoving observers that Fermi- 
propagate their spatial reference frames) of disturbances in the matter and the gravitational field, 
represented as wavefronts by the characteristic 3-surfaces of the system, are obtained. In particular, 
the Weyl curvature is found to account for two (non-Lorentz-invariant) Coulomb-like characteristic 
eigenfields propagating with v = and four transverse characteristic eigenfields propagating with 
\v\ = 1, which are well known, and four (non-Lorentz-invariant) longitudinal characteristic eigenfields 
propagating with |w| = |. The implications of this result are discussed in some detail and a parallel 
is drawn to the propagation of irregularities in the matter distribution. In a worked example, we 
specialise the equations to cosmological models in locally rotationally symmetric class II and include 
| the constraints into the set of causally propagating dynamical variables. 
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The natural way to relate gravitational fields to physical measurements is by detecting their tidal and their volume- 
changing effects: monitoring the relative accelerations between a set of test particles in free fall, measuring the 
gravitationally induced deformations between the ends of a semi-rigid detector, or detecting the distortions of the cross 
section of an infinitesimal beam of light that travels freely through spacetime. This viewpoint has been emphasised 
strongly by a number of authors, particularly, Pirani ]4l[] , Szekeres [Q, and Wheeler [p3| . In mathematical terms, 
the connection between a gravitational field and the measurement process is made through the geodesic deviation 



X 



equation ('GDE') , introduced into general relativity by Synge 47 1 , wherein the Riemann curvature of a given spacetime 



5^ , manifold ( M. , g ) provides the force term responsible for the worldlines of test matter to deviate from the straight 
&)[)■ lines of Minkowskian geometry.]] If the (timelike or null) tangents to these worldlines are denoted by V a , the second 
rate of change along V a of the (orthogonal) connecting vectors rj a is determined by the Riemann curvature ( see, e.g., 
Ref. |H): 



V c V c (V b V bV a ) = - R a bcd V b if V d . (1.1) 

Geometrically the Riemann curvature is often defined as the linear operator which encodes the non-commutativity of 
two covariant derivatives on (M, g), i.e. 

2V [a V b] v c :=R ab c d v d , (1.2) 

for any vector v a , from which the above equation follows. As is standard, the Riemann curvature tensor can be 
decomposed into the Weyl curvature tensor, the Ricci curvature tensor, and the Ricci curvature scalar according to 

Rabcd — Cabcd + R a [c 9d]b ~ Rb[c9d]a ~ \ R9a[c9d]b ■ (1-3) 



*Electronic address: henk@gmunu.mth.uct.ac.za 
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: Or through the closely associated relative deviation equation for particles subject to non-gravitational forces, and, hence, 
moving non-geodesically (see, e.g., Ehlers ); however, the essential gravitational effects are already incorporated in the 
GDE, so we concentrate on that equation in this paper. 



1 



Consequently these are the force terms in the GDE. If, as is the case in relativistic cosmology, there exists an invariantly 
defined future-directed unit timelike vector field u, the Weyl and the Ricci curvatures of (M, g, u) can be further 
decomposed with respect to the group of spatial rotations in a 1 + 3 covariant way |ll||l5| according to 

C ab cd (u) =4u^u [c E^ d] +4h^ [c E b ^ d] + 2e abe u [c H d]e + 2e cde u^H b ^ (1.4) 

-Rab(u) = \ (p. + 3p - 2A) u a u b + 2 q (a u b ) + \ (p, - p + 2A) h ab + ir ab , (1.5) 
respectively, where in Eq. ( |1.4| ) E ab and H ab are the electric and magnetic Weyl curvature tensors (relative to u), 



has triggered a growing 



and in Eq. (L5) the Einstein field equations ('EFE') R ab — i Rg a b + ^-9ab = T ab have been used to equate the Ricci 
curvature to the energy-momentum-stress of the matter source fields, T ab (u) := fi u a u b + 2 q^ a u b ^ + p h ab + n ab . The 
matter fluid variables /i, q a , p, and ir ab are the energy density, energy current density, isotropic pressure, and tracefree 
anisotropic pressure (relative to u); h ab :— g ab + u a u b is the projection tensor into the instantaneous rest 3-spaces 
orthogonal to u with 3- volume element e abc := — e de f g h d a h% h' c u 9 . 

It is clear, then, that all these terms enter into the GDE and can have measurable physical effects; hence, their 
propagation behaviours are of physical importance. Tying the 1 + 3 viewpoint to an orthonormal frame description 
which includes the full Riemann curvature, the dynamical equations for the associated geometrical field variables are 
provided by the frame basis commutation relations, the Jacobi, the Ricci and also the Bianchi field equations | pl| . 
These determine the propagation modes of all the above quantities, but not immediately in a mathematically most 
desirable way, because the hyperbolic nature of the overall equation system is not cle arly ap parent. 

Recently research activity in the arena of numerical relativity ( see, e.g., Refs. |^,[l|, ^2| , p2||43 
interest in formulating evolution systems for dynamical field variables in such a way that their propagation velocities 
transparently respect relativistic causality, and where physically motivated boundary conditions can be implemented 
more easily than in the standard ADM 3 + 1 formalism Q . A natural framework for this purpose is provided by the 
first-order symmetric hyperbolic ('FOSH') form for sets of (quasi-linear) partial differential equations, that Friedrichs 
introduced in the mid- fifties [E6[ . Fortunately, most sets of dynamical equations occurring in theoretical physics can 
be cast into such a form p7[ff For a physicist this is especially attractive for at least two reasons: first, there do 
exist mathematical proofs that show (local) existence and uniqueness of solutions to the related Cauchy initial value 
problem; secondly, FOSH evolution systems allow for the propagation of data that are non-analytic across certain 
(characteristic) 3-surfaces, which can be interpreted as wavefronts. A standard reference on both aspects is the book 
by Courant and Hilbert p|. 

Suppose a collection of dependent field variables u A — u A {x^), which are functions of a set of local spacetime 
coordinates {x M }, satisfies the equation system 

M AB ^d^u B = N A . (1.6) 

The objects M AB ^ = M AB ^{x v , u c ) and N A = N A (x^ ,u B ) denote four symmetric matrices and a vector, respec- 
tively, each acting in a sp ace of dimension equal to the number of dependent fields. Conventionally one refers to the 
left-hand side of Eq. ([lj) — the 'principle part' — as representing dynamical interactions between certain physical 
fields, the right-hand side representing kinematical ones (cf. Geroch ). The set of equations ( p"lj| ) is hyperbolic if 
the contraction M AB with the coordinate components of some past-directed timelike 1-form n a yields a positive- 
definite matrix; it is causal if this contraction is positive-definite for all past-directed timelike 1-forms n a pl[ ]. If it 
satisfies all these conditions, it is said to be a FOSH evolution system. 

In the present work we inquire into the optimal formulation of a FOSH evolution system for cosmological models 
(M, g, u) with barotropic perfect fluid matter source fields, starting from the dynamical equations of the extended 
1 + 3 orthonormal frame ('ONF') formalism which includes the Weyl curvature variables, as discussed recently in Ref. 
pl| . The fact that the Bianchi field equations for the Weyl curvature, which are of third order in the derivatives of the 
metric g, naturally lead to FOSH equations, was already recognised by Friedrich in the early eighties p3| . Moreover, 
lately he employed them in work on the derivation of FOSH evolution systems for vacuum gravitational fields |24|] 
and for gravitating perfect fluid bodies [^5| . Our choice of dependent field variables as provided by the extended 1 + 3 
ONF formalism shares the same philosophy as Friedrich's, but additionally we have fully decomposed the spacetime 
connection into its irreducible parts (hence, obtaining more transparent equations), and we introduce a set of local 



2 Incidentally, Geroch leaves the case of the (vacuum) EFE in his nice work open; as other papers cited in this introduction 
and the present work show, a FOSH evolution system for general relativity can be obtained, e.g., by either employing the Bianchi 
field equations or introducing auxiliary variables based on first derivatives of the metric. 



2 



coordinates by following the 1 + 3 threading picture |3^|| • The latter is of particular relevance when the timelike 
reference congruence of the setting has non-zero vorticity. 

Special attention in our investigation will be given to the properties of a longitudinal component in the Weyl 
curvature as described by Szekeres in work on the propagation of gravitational fields in matter ||48| , |49f . Friedrich (in 
the vacuum case) discounted the propagation modes for these components as not being physically significant, because 
the characteristic 3-surfaces along which they travel are not Lorentz-invariant j24J]. We suggest an opposite view: 
that despite this feature, these modes must be seen as physically significant because of their detectable effects as 
determined by the GDE. We briefly consider some consequences of this suggestion. 



We specify our conventions, choice of units, and some of our notation in appendix A 1 



II. GENERAL 1 + 3-DECOMPOSED DYNAMICAL EQUATIONS 

For easy reference, we dedicate this section to a complete exposition of the 1 + 3 ONF dynamical equations, where 
the energy-momcntum-strcss of the matter source fields is assumed to have the form of a general dissipative relativistic 
fluid. As described in detail in Rcf. |2l|, for a cosmological model ( M, g, u) a set of four orthogonal unit basis fields 
{ e a } is chosen and its timelike member eo is identified with the unit tangent u to the matter flow lines (here as 
viewed from the 'particle frame' perspective). This fixing implies that we reduce the initial six-parameter Lorentz 
transformation freedom in the choice of { e Q } to a three-parameter freedom that allows for arbitrary rotations of the 
spatial frame { e a } only. 

A. Connection components and commutators 

The 24 algebraically independent frame components of the spacetime connection T a i, c can be split into the set^| 

TaOO = U a = TfoOo (2-1) 
r a 0/3 = 3 © 0~af3 + 0~ a p — tap-y — Tp a /3 a (2.2) 

r a /30 = £af3y ^ = 0af3 (2.3) 

r Q/ 3 7 = 2 ci[ a (5 (3 ] 7 + e 1 s[a nS p] + \ tups n 5 j = T F p ia . (2.4) 

It contains the familiar kinematical fluid variables: u a , its relativistic acceleration, 0, its volume rate of expansion, 
o~ a p — o-p a , its rate of shear (with a a a — 0), and cu a , its vorticity. Moreover, Q a is the rate of rotation of the 
spatial frame { e Q } with respect to a Fermi-propagated basis. Finally, a a and n a p = ng a are 9 spatial commutation 
functions. 



With our definition of the covariant derivative on (Ai, g, u), given in appendix Al, and Eqs. (2.1) - fl2.4|) , the 
condition that the spacetime connection be torsion-free, Vr a Vw/ = for any spacetime scalar /, translates into the 
commutator equations 

[eo, e Q ] (/) - u a e (/) - [ | 6 6 J + a J + e/ 7 (^ + W) } e p (f) (2.5) 
[e a ,e^] (/) = 2e Q/ 3 7 w 7 e (/) + [ 2a [a 5 1 fi] + e aj 3&n &1 ] e 7 (/) . (2.6) 



B. Constraints 



The following set of relations does not contain any frame derivatives with respect to eo. Hence, it is commonplace 
to refer to these relations as 'constraints'^] From the Ricci identity for u as well as the Jacobi identity we have the 



3 To facilitate comparison, we indicate how these frame components of the connection relate to the respective variables used 
by Friedrich in Ref. J5^|, which we denote by IV abc- 

4 Even though this terminology is problematic in the generic case when eo = u has non-zero vorticity and local coordinates 
are introduced according to Eq. ([3.6|) below. 
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D;,cr ofc -equation (Ci) Q , which, in Hamiltonian treatments of the EFE, is also referred to as the 'momentum constraint', 
the D a w a -equation (C2), and the ff a h-constraint (C^)"' 3 , respectively; the once-contracted second Bianchi identity 



yields the DbE ab - and Dfci/ ab -equations (d)" and (C$) a J_5|Jl7|; the constraint (Cj) a again arises from the Jacobi 
identity, while, finally, (Cc) a ^ and (Cq) stem from the EFE. In detail, 

0=(d) a := {ep-iap){a aP )-l5 a0 e p {Q)-n a p J j + q a 

+ e a ^ [(e + 2up-ap) - n ps a% ] (2.7) 

= (C 2 ) := (e Q -u a -2a a )(w a ) (2.8) 

0=(C 3 ) Q/3 := i? Q/3 + (<5 7<Q e 7 + 2u (Q + a <Q )(w /3> )-in\a Q/3 + 3n< Q 7( j /3 >T 

- e 7d > [ (e 7 - a 7 ) (<A) + n% co 5 } (2.9) 

0=(C 4 ) tt := (ep-3ap)(E a ' 3 + ^n a ' 3 )-^6 a ' 3 e p ( f i) + leq a -^a a pq^-3u;pH^ 

- e a ^ [ apt H% -§ W/3 « 7 + n/sa (£\ + § ^ 7 ) ] (2.10) 

0=(C 5 ) Q := (e^-3a^)( J ff Q/3 ) + ( M + p)^ + 3^(^-i^)-in^g' 3 

+ e Q ^ [ § ( e/3 - a/3 ) (g 7 ) + ^5 (£ 5 7 + i ^ 5 7 ) - n/M # 5 7 ] (2.11) 

= (Cj) a := (ep-2 ap) (n al3 ) + § 6 uj a + 2 a a p lj 13 + e a ^ [e (a 7 )-2 ojp ft 7 ] (2.12) 

0=(C G ) Q/3 := *S aP + ±&<r al3 -a^a^ -J a uj® + 2J a fL® ^ (E al3 + \ir al3 ) (2.13) 

= (C G ) := *i?+§e 2 -(aV^ + 2Kw Q )-4(w Q f! Q )-2 M -2A , (2.14) 

where 

*S Q/ 3 := e (Q (a /3 )) + 6( Q/3 ) - e 7 " 5 <Q (e| 7 | - 2o| 7 |) (ng) 5 ) (2.15) 

*R:=2(2e 0t -3a«)(a a )-§& O a (2.16) 

& Q/3 := 2 ?i Q7 n 7 ^ - n l 1 n a p , (2-17) 

and angle brackets denote the symmetric tracefree part. If oj a = 0, such that the matter flow tangents u become 
the normals to a family of spacelike 3-surfaces T: {t — const}, the last two constraints in the set correspond to the 
symmetric tracefree and trace parts of the once-contracted Gaufi embedding equation. In this case, one also speaks 
of (Cq) as the generalised Friedmann equation, alias the 'Hamiltonian constrain t' or the 'energy constraint'. When 
cu a = 0, the set of constraints should also include the commutation relation Eq. fl2.6|). 



C. Evolution of spatial commutation functions 

The 9 spatial commutation functions a a and n a p are generally evolved by Eqs. (40) and (41) given in Ref. @; these 
originate from the Jacobi identity. Employing each of the constraints (Ci) a to (C^)"' 3 listed in the previous paragraph, 
we can eliminate e a frame derivatives of the kinematical fluid variables 0, a a p and uj a from their right-hand sides. 
Thus, we obtain the following equations for the evolution of the spatial commutation functions: 

e (a a ) =-i(6^-| G a p) (vP + a?) + \ n a p ^-\q a 

~ \ e a ^ [ (up + ap) - np S a% -(ep+up-2 ap) (fi 7 ) ] + \ (C x ) a (2.18) 

e (n al3 ) = - \ 6 n af3 - ct (q 7 n®"* + \ a af} n 7 7 - (u<" + a {a ) oj® - H al3 + (5~< {a e 7 + u <a ) (Q® ) 
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%8^[2 (ii 7 + a 7 ) - a~> s n% + (e 7 + u 7 ) (ST) ] 

e^< Q [ (u 7 + a 7 ) <A - K + 2 fi 7 ) rA ] - | (C 2 ) + (Ca)"' 3 . (2.19) 



D. Evolution of kinematical fluid variables 

The evolution equations for the 9 kinematical fluid variables O, o~ a p and uj a are provided by the familiar Ricci field 
equations, i.e. 

e (9) - e a (u a ) = -\<d 2 + (ii a -2 a a ) ii a - (a a pa a ) + 2 (u a uj a ) - \ (fi + 3p) + A (2.20) 

e (a al3 ) - <F< a e 7 (u /3) ) = - § 6 a afi + {u {a + a (o ) u"> - cr< a 7 - Lu<~ a u 0) - (E afs - \ n a(3 ) 

+ e~* s{a [ 2 n 7 o® s - vP\ us } (2.21) 



(u a ) - | e Q/57 ep(u 7 ) = - | 6 cu a + u a p uj p - \ n a p u 13 - \ e a ^ [ ap u 7 - 2 Sip u 7 ] . (2.22) 



e ( 



E. Evolution of Weyl curvature and matter variables 

Finally, we have the Bianchi field equations for the 10 Weyl curvature variables E a p and H a p and the 4 mat- 
ter variables p and q a , which are obtained from the once-contracted and twice-contracted second Bianchi identity, 
respectively p|P,p5| : 



e (E afl + \ ir a0 ) - e^ {a e 7 (H® s ) 

+ \ <F <Q e 7 (/> ) = - § (fi + p) a a P - 9 [E a!i + § 7r a/3 ) + 3 (E^ - ± tt^) 
+ § n\ H a(i - 3 n% H®~< - ± (2 + a {a ) q (3) 
+ e^ {a [ (2u 7 -a 7 )H® s 

+ K + 2 fy) + | tA) + i n« 7 % ] (2.23) 

e (H al3 ) + e^ a e 7 (E® s - \ n f) \) = -OH afj + 3 a< a 7 + § J a q® 

- \ n t 1 {E al3 - \ Ti af) ) + 3 n% (E^ - \ tt^t) 

+i^> 7 g 5 + (cj 7 +20 7 )ff^ ] (2.24) 

e (<?<*) + 5^ e/9 (p) + e^(Tr^) = -| 6 g a - / - (/x + p) u a - (up -Sap) tt^ 

- e Q/ ^ [ (up - Sip) q 7 - np 5 n 5 7 ] (2.25) 

e (n)+e a (q a ) = -Q( f i+p)-2(u a -a a )q a - (a a p^ a ) . (2.26) 



The general extended 1+3 ONF dynamical equations (2J5) and ( 2.18 ) - ( |2.26j ) do not directly form a FOSH 
evolution system, for two reasons. First, they do not provide evolution equations for any of the geometrical variables 
m q , Sl a , p and 7r Q ^J^] In the case of Sl a this is a reflection of the freedom of choice of a particular frame { e a }. The 



5 The time derivative of the latter is indeed included in one of the Bianchi field equations, but not in a 'pure' form; because 
of the dynamical meaning of 7r a( g discussed below, we regard this as an equation for E a/ 3. 
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indeterminacy of propagation of the remaining three variables reveals the necessity of tying the description of non- 
vacuum gravitational phenomena to a thermodynamical description of the matter source fields. This may be achieved, 
for example, by a phenomenological scheme for dissipative relativistic fluids or a relativistic kinetic theory approach. 
Second, the combinations of derivatives occuring in the above equations do not have the required symmetric structure. 
Given our present goal, the target is to choose a suitable matter description and then shuffle the evolution equations 
for the full set of variables into a FOSH form. 



F. New set of frame variables for rank 2 symmetric tracefree tensors 



In practical applications of the 1 + 3 ONF formalism it proves helpful to introduce a new set of variables adapted 
to the tracefree condition A a a = and the mag nitude A 2 := \ A a l3 A' 3 a > of any 1 + 3 invariantly defined rank 2 
symmetric tracefree tensor field A orthogonal to u. In our case we have, in particular, A a p € { a a p, E a p, H a p, Tr a p }. 
This new set of irreducible frame components is defined by ( cf. Ref . [El] ,|21 ] ) 



A+ :=-§ A n = § (A 22 



At := V3A 



•2:i 



A 33 ) A_ := & (A 22 - A 33 ) 
A 2 := V3A 31 



-4, 



(2.27) 



It has the property that the magnitude of A assumes the explicit form 

A 2 = i [ {A+f + (A_) 2 + (A,) 2 + (A 2 ) 2 + (A3) 2 ] ■ (2.28) 

It should be emphasised that this decomposition arbitrarily adapts to the spatial ei-axis. However, this is only a 
matter of convention, and by a cyclic permutation of indices 1 — > 2 — * 3 — > 1 one can easily adapt to any of the other 
spatial axes. 

In principle, one can extend this procedure to also include the spatial commutation functions n a p, i.e., splitting 
them into a trace and a tracefree part, where the latter is further subdivided according to the scheme above. However, 
as n a p has no immediate invariant meaning in the 1 + 3 covariant picture, this possibility will not be pursued any 
further here. 



III. (M, G, U) WITH BAROTROPIC PERFECT FLUID MATTER SOURCE FIELDS 



We now specialise the energy-momentum-stress tensor of the matter source fields appearing in the dynamical 
equations of section || to a perfect fluid form (with matter flow tangents u) by setting 

= q a =n a p (3.1) 

throughout,^ and we assume a barotropic equation of state, i.e. , a functional dependence of the isotropic fluid pressure 
on the energy density only, 

p=p(p). (3.2) 

We define c 2 (/i) : = dp{p)/dfi as the isentropic speed of sound. Under the assumptions ( |3.l| ), the evolution equation 
for q a , Eq. ( 2.25 ), reduces to the new constraint 

= (C PF ) a :=5 al3 e (p) + (fi+p)u a , (3.3) 

which is often called the momentum conservation equation. 



We also discard the cosmological constant, A = 0. 



G 



A. Derivation of a FOSH evolution system 



By applying the commutator (2.5) to / = p and using Eqs. (|3.2|), (p\3|), ( |2.26 ) and (2.7), we derive an evolution 
equation for the acceleration u a of the matter fluid tangents u given by 



e {u a ) - c 2 s b p {\ 6 6 aP + a al3 + e a ^ w 7 ) 



d 2 p 



(p+p)-c 2 s + ±}eu a - (up + Zc 2 s ap) o ali - c 2 s n a p u p 



_l_ e "/37 [ (2c 2 — 1) uip u 7 + flpii-y — c 2 ap lo 1 — c 2 np$ a\ 



(3.4) 



Here we have suppressed an evolution equation for (Cpf) u whic h ar is es alo ng th e wa y. Next, contracting the com- 
mutator (2.6), again applied to / = p, with e a/31 and using Eqs. (3.2), (2.26) and (3.2) leads to the identity^ 



= e a ^ (ep - a ) (u 7 ) - n a p /-2^9u a 



(ji + p)- 1 [ e a ^ (ep - ap) (CW) 7 - n a p (C PF f + (c7 2 + 1) e a ^ up (C PF \ 



(3.5) 



This identity constitutes the key step in achieving FOSH form for the evolution subsystem ( p. 20 ) - (2.22) and (3.4) 
that links the kinematical fluid variables 0, a a p and u> a to u a and establishes the sound cone structure on ( M., g, u). 
The trick is to add, on using Eq. (3.5), e Q ^ 7 ep(u 7 ) to the left-hand side of Eq. ( 2.22| ), i.e., to change its principle 
part to the new form eo(iv a ) + h € a '^~ep(u 1 ). 

Having adopted the Lagrangean viewpoint for the description of the dynamics by making the identification eo = u, 
we now proceed to introduce a set of (dimensionless) comoving local coordinates { x^ } = {t, x l }. The natural way 
of doing this is via the 1 + 3 threading approach discussed in detail by Jantzen et al |32| and Boersma and Dray Jm. 
Defining coordinate components of the e a by e Q M :— e a (x tJ '), the frame spanning basis fields can be expressed as [£lj 



e = eo' 1 d„ := A/" 1 d t 



= e a »d„ :=e a l (Af< fit + 



(3.6) 



where M = M(t,x l ) is the threading lapse function and Midx 1 = Mi(t,x 3 ) dx l is the dimensionless threading shift 
1-form. The inverse of the threading metric is h lJ :— 5 a ^ e a l ep 3 . Note the presence of the temporal partial d eriv ative 
i n th e 'spatial' frame fields e a . With the expanded frame fields ( p\6| ) inserted, the commutator equations (2.5) and 
(|1| now yield @ 



e a l [ d t Mi + M- 1 (Mi dtM + d t M) } = u 



(3.7) 



M- 1 d t e a l = - [ \ 6 5 J + a J + ej 7 + W) } 



(3.8) 



M e [a l e f (Mi d t M 3 + diM 3 ) = e aPl u/> 



(3.9) 



2 e [a l { M t dtepf + d^f ] e 7 , = 2 a [a 6~>p ] + e a p S n* 7 



(3.10) 



where e a i is defined through the relation e a ie a 3 = 5 3 i. The threading shift 1-form will be non-zero, Mi ^ 0, if the 
worldlines of the matter fluid are rotating, i.e., u> a 0. In this case the local rest 3-spaces do not mesh togeth er to 
form a fa mily of spacelike 3-surfaces T: {t — const} everywhere orthogonal to u. Substitution of dt Mj in Eq. (3.9) 
from Eq. (3.7) generates a further constraint relation to be satisfied. 



One way of picking the threading lapse function M is provided by parameterising the integral curves of u by physical 
proper time. Then we have 



7 In 1 + 3 covariant terms this is just = e abc D b u c - 2 c 2 3 6 cu a - . . . (see Ref. @ , Ch.2 ) 
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e ' 



U " = M - 1 ^o , 



(3.11) 



where Mo denotes a constant threading lapse function of unit length. If, furthermore, the spatial frame { e a } is 
chosen to be Fermi-propagated along u, i.e., 



n a = , 



we obtain from Eqs. (3.7) and ( |3.8| ) reduced evolution equations for Mi and e a l ; 

<9 t Mi = u a e a i = Ui 



(3.12) 
(3.13) 



M^dteJ = -[ hQSj+aJ 



(3.14) 



An alternative to a constant threading lapse function and a dynamical threading shift 1-form is the prescription for 
barotropic perfect fluids given by Salzman and Taub W4|, where 



M = M exp 



dp' 



(m + pO 



M i = M i {x j ) ; 



(3.15) 



here M is dynamical while Mi is coordinate time independent. This choice reduces the number of dependent field 

F""" 1 "' 1 " tl,E — " """""" ™ '' v ! "" : " 



Our preparations for the derivation of an evolution system in FOSH form for barotropic perfect fluids from the 
general 1 + 3 ONF equations are now complete. To this end we will use the set given by the matter content specifying 
equations ([O) and (|3.2| ), the frame fixing equa tions (3.11) and ( 3.12| ), the evolution equations ( |3.13 ), (3.14), ( 2.18) ) 
- ( [2.24 ), ( |2.26| ) and ( |3.4[ ), and the identity (3.5). Of imm edia te interest for the FOSH structure is only the principle 
part of the evolution equations, the left-hand side in Eq. (1.6), that describes the dynamical interactions between the 
various fields. In terms of the frame derivatives e a we can represent it by 



The frame derivative principle part: 



M A " a e a (u B ) 



(3.16) 



Using tracefree- adapted irreducible frame variables as defined in Eq. ( 2.27 ) for each of the fluid rate of shear and 
the electric and magnetic Weyl curvature, a FOSH evolution system can now be obtained, taking certain linear 
combinations of the equations where necessary, for the following set of 44 dependent dynamical fields: 



The dependent geometrical field variables: 



( H frame ^ 

^3-con 

"few,! 

U>kin,2 

"few, 3 

Urnat 
\ UWeyl / 



u> frame 
^3-con 

Wfe m ,2 

^mat 
UWeyl 



e-cx. i 

a a , n al3 



Mi) T 

T 



tlx, |(© - 2(7+), M-£T 3 + W 3 ), (775^2 - UJ 2 ) 



[ U 2 , |(9 + <J + + \/3(7_), (^CTl + Wl), (77503 - UJ 3 ) 



V3 L 



(3.17) 



= [ u 3 , |(0 + cr + - V3t7-), (775^2 + lo 2 ), (77501 - w x ) ] 



V3 L 



[ E + , £L, Ex, E%, E 3 , H + , Hi, H 2 , H 3 ] J 



The symmetric (44 x 44)-matrices M ABa occurring in Eq. (3.16) are found to assume the forms 
The matrices: 

\ /o 12 



/ ll2 



M 



ABO 



K 



K 



K 



M 



AB 1 



o 9 



B; 



Bo 



(3.18) 



CiJ 



8 



/Oi 



M 



AB2 



B, 



M 



-jAB 3 



B 2 



C 2 J 



B 2 



B a 



Bi 



(3.19) 



C 3 J 



K := 



B, := 



Bo 



Vl 



, Ba :— i 



(3.20) 



C a := 



-1/2 



1/2 



1/2 



1/2 



(3.21) 



V 



V3/2 -1/2 



V^/2 1/2 



V3/2 . \ 
-1/2 . 

1/2 



V3/2 
1/2 



-1/2 



1/2 



1/2 



(3.22) 



C, := 



\-V3/2 -1/2 



-1/2 



1/2 



V3/2 
1/2 



V3/2 1/2 



1/2 



-1/2 



-\/3/2\ 
-1/2 



(3.23) 



With the exception of the submatrices K and to B 3 in the kinematical branch of our FOSH evolution system 
that contain c 2 s = Cj(/i), the remaining components of the M ABa are just numerical constants, since we evaluated the 
principle part in the 1 + 3 ONF b asis { e a } (and because of reasons which we will briefly come back to below) . When 
the frame derivatives e a in Eq. (3^16) ar e sub stituted for in terms of their coordinate components and the partial 
deriv ative s introduced in Eqs. (3^) and ( |3.1l| ), the non-constant components of the symmetric matrices M AB ^ in 



Eq. (1.6) with respect to the local coordinate basis can be easily read off. Finally, hyperbolicity will follow 

from the reality of the eigenvalues we obtain when determining the characteristic 3-surfaces in the next subsection. 







B. Characteristic 3-surfaces, propagation velocities and eigenfields 



The set of characteristic 3-surfaces {</> = const} underlying a FOSH evolution system can be interpreted as a 
collection of wavefronts with phase function cj> across which certain physical quantities may be discontinuous. The 
associated characteristic eigenfields propagate along so-called bicharacteristic rays with tangents k within these 3- 
surfaces at velocities v, where v — tanha is a function of the hyperbolic angle of tilt a between k and u fl2| ]. To 
determine the characteristic 3-surfaces, we can assume without loss of generality that at a particular spacetime event 
the 1 + 3 ONF is oriented in such a way that the orthogonal gradient vector fields £ a := V a </>, are tangent to a 
2-surfaces spanned by u and ei. This choice of frame, which then only leaves the freedom of spatial rotations about 
the ei-direction, can be made due to the local isotropy of the characteristic cones (see, e.g., Refs. §1311)). Hence, 
locally the frame axes and the given coordinate lines are made to coincide. In particular, at a point we can take 
coordinate components of ei where = ei 2 = ei 3 and also Mi = 0;|] however, generically their time derivatives will 
be non-zero. Thus, we obtain 

e Q = M^d t , e 1 =e 1 1 8 1 . (3.24) 



1. The characteristic condition 

The characteristic condition the vector fields £ have to satisfy is 

= det[M AB "£ M ] . (3.25) 

With the choice of frame outlined, their coordinate components are given by £ M = —vu^ + e 1 ^ fl^ , |29| , |3"lf , where 
the parameter v coincides with the different possible propagati on v elocities of the characteristic eigenfields^ and e 1 ^ 



denotes the inverse coordinate components of ei. With Eq. ( 3.24 ) we find, therefore, £ M = vMqS ^ + (ei 1 ) 1 # M , 
leading to 

= dot [vM M AB + (ei 1 )" 1 M AB 1 ] 

= (M ) 44 4 v 30 (v - c s ) 3 (v + c s ) 3 (v - l) 2 (v + l) 2 (v - \f (v + i) 2 . (3.26) 

Clearly all roots v of this equation are real (establishing hyperbolicity; cf. Ref. |Q ), and no root v is single- valued, 
implying a set of characteristic 3-surfaces each of which is degenerate. 

2. Characteristic velocities and eigenfields 

First, we find 30 characteristic eigenfields of our FOSH evolution system that propagate with velocity v\ = with 
respect to observers comoving with u; these are 

uf x) = [e a \ M h a a , n a0 , | (9 + a + ± ^<r_ ) , (^-a 1 ±uj l ), (^<t 2 -lj 2 ), (-^<j 3 +u; 3 ), ft, E+, H+ } T . (3.27) 



Second, upon diagonalisation of the ei-adapted, reduced principle part in Eq. (F6), the set of 6 characteristic 
eigenfields propagating with velocities «2,3 = ± c s along the sound cone is found to be 

u f 2 ) = [ -^( c su 1 ±l(e-2a + )), -jl==(c s u 2 ±(±o- 3 -u> 3 )), -^=( Cs u 3 ±(^a 2 +u J2 ))} T . (3.28) 
Here, and in the following, the upper sign applies to outg oing modes and the lower one to incoming modes.^ By 



construction, the tracefree-adapted variables of Eq. (2.27) clearly exhibit the purely longitudinal character of the 



The vorticity of the matter flow lines cannot be felt at a single point. 
9 It also gives the slope of £ with respect to ei. 

10 In observational cosmology, where we are faced with a situation of 'here and now' in spacetime, it is the incoming modes 
that are of physical relevance. 
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first two eigenfields)^] and the semi-longitudinal of the latter four with respect to the (assumed) spatial propagation 
direction ei. In the pressure-free case, p — =>• u a — 0, c s = 0, the remaining kinematical fluid variables 0, a a f3 and 
uj a also propagate with v\ = Q.p| 

Third, the set of 4 characteristic eigenficlds propagating with velocities U4 5 = ± 1 along the null cone is 



l (3) 



[^j(^-=FHi),^(^±H-)P 



(3.29) 



each pair corresponds to one of the two poss ible polarisation states of the freely propagating gravitational field. 
Again, the tracefree-adapted variables of Eq. (2.27) nicely reveal that these eigenfields are purely transverse to the 
(ass umed) spati al pr opagation direction ei. The full propagation equations of these modes are given in appendix 
A 2 ; in appendix A 4 one can find expansions of all Weyl curvature characteristic eigenfields using the gauge-invariant 
variables of Bardeen for linearised perturbations of FLRW cosmological models [|| . 

Finally one obtains 4 characteristic eigenfields propagating with velocities U6.7 = ± \ along timelike 3-surfaces 
which are 



>4) = [^(E 3 TH 2 ), -L(E 2 ±H 3 )]" j 



(3.30) 



The tracefree-adapted variables of Eq. ( 2.27 ) classify them as longitudinal modes of the Weyl curvature, which 
also come in two possible polarisation states. They have be en i dentified before by, e.g., Szekeres p8| , p9[ . The full 
propagation equations of these modes are given in appendix A 2. Their physical relevance will be further discussed 
below. 



The results of Eqs. (3.27) - ( 3.30| ) lead to a number of observations 
local causal structure of (Ai, g, u) (and has coordinate components constructed by 



The spacetime metric g, which embodies the 
77a&e a AI e fc „), propagates 
= ( see also Ref . E2| ) . Parts of the spacetime connection 
0, while the remaining parts follow the sound cone. With 



along the timelike reference congruence u itself, i.e., with v 
(which contains first derivatives of g) also propagate at v = 
the present geometrical set of dependent field variables, it is only within the Weyl curvature (which is of second order 
in the derivatives of g) that one finds modes that propagate changes in the state of a gravitational field at the speed 
of light. 

As to th e ma gnitud es of the different propagation velocities it should be noted that, by use of the Bianchi field 
equations ( 2.23 ) and ( 2.24 ), causal propagation of the gravitational field modes with \v\ £ { 0, |, 1} falls out auto- 
matically; no further assumptions are required. To ensure causal propagation of pressure perturbations in the matter 
fluid, on the other hand, we need to impose the condition < c s < 1. 

For vacuum gravitational fields as well as perfect fluid spacetime geometries (in matter comoving description) , the 
different characteristic 3-surfaces associated with the Weyl curvature do not intersect due to the 'linear degeneracy' 
of the principle part in the Bianchi field equations. This means it is not genuinely non-linear and, therefore, discon- 
tinuities ('shocks') cannot form in the process of evolution. The only Weyl curvature discontinuities that propagate 
are those contained in the initial data. The situation is different in the sound cone branch of our FOSH evolution 
system, where the (non- const ant, equation of state dependent) characteristic velocity c s explicitly occurs in the related 
eigenfields given in Eq. ( 3.28 ), and shock formation is possible ( on these issues see, e.g., Bona in Ref. §). A different 
phenomenon altogether are the 'shell-crossing singularities'. At these events the Lagrangean description breaks down 
because integral curves of u intersect. 



3. Polarisation states and the NP formalism 

We mentioned before that by our method of choosing a 1 + 3 ONF {e a }, with e aligned along the uniquely 
defined matter flow lines and ei locally aligned along an arbitrarily chosen propagation direction (for the purpose 
of determining the characteristic structure of our FOSH evolution system), the remaining frame freedom, given this 
choice of propagation direction, consists of spatial rotations about the e!-direction through an angle ip only. Focusing 
on the characteristic eigenfields of the Weyl curvature, this implies that out of, say, the 4 outgoing modes, in general 



11 The quantity |(6 - 
12 When u a = 0, Eqs. 



2cr+) is just the longitudinal component |0 + 011 of the rate of expansion tensor. 

(2.20) - (2.21) do not need to be multiplied by c\ to obtain a symmetric structure for the principle part 



in the kinematical branch. 
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we can set only one to zero, as the following transformation equations clarify: 
Transverse v — ±1 modes: 

(E- + Hi) = (£L + Hi) (cos 2 ip - sin 2 ip) - 2 (E x ± HJ) cosip sin 93 (3.31) 
(Ei ± #_) = (Ei ± (cos 2 93 - sin 2 93) + 2 (E- + #1) cos^ sin ip ; (3.32) 

hence, the magnitude of the difference between the two angles that set one or the other polarisation state of the 
transverse component to zero is given by \(fi — p 2 \ = \. 

Longitudinal v = ± 5 modes: 

(E 3 + H 2 ) = (E_ + Hi) cos 95 - (E 2 ± H 3 ) simp (3.33) 
(E 2 ±H 3 ) = (E 2 ±H 3 ) cos ip + (E 3 + ff 2 ) sin y> ; (3.34) 

in this case, the magnitude of the difference between the two angles that set one or the other polarisation state of the 
longitudinal component to zero is now given by \pi — p 2 \ — ■§■ 

Concluding the current considerations on the features of the characteristic eigenfields of the Weyl curvature, for 
comparison we can construct the complex scalars introduced by Newman and Penrose |4(J] in the null tetrad canonically 
associated with our 1 + 3 ONF. With the help of Eqs. (C.26) and (C.27) given in Ref. u%, we find (in the language 
of Szekeres [ ^8p9[| ) for the incoming and outgoing transverse radiative parts, 

V =-^(E_+Hi)-i-^(Ei-H_) and tf 4 = (E_ - H t ) + i ^ (E x + H_) , (3.35) 
for the Coulomb part 

* 2 = - \E + - i § H+ , (3.36) 
and for the incoming and outgoing longitudinal radiative parts 

^^ = -^^ + H ^ +i ^T^ E ^ H ^ and VE '3 = ^(^3-^2)+^(i?2 + ^3) ■ (3-37) 



C. Propagation of constraints 



The analysis of the propagation of the constraints (2/7) - ( 2.14 ) is a slightly delicate issue, as in the given context the 



tangents u to a generic matter flow have non-zero vorticity; hence, there exist no well-defined spacelike 3-surfaces T: 
{t = con st) everywhere orthogonal to u, and the frame derivatives e a contain a temporal partial derivative according 
to Eq. (|3.6|) . In principle this problem could be circumvented by slicing (Ai, g, u) with any arbitrary family of 
spacelike 3-surfaces T: {t = const}, on which one alternatively evaluates the constraints and propagates them along 
the timelike normals (which will be tilted with respect to u) . In this non-comoving approach the energy-momentum- 
stress tensor would appear imperfect. A preliminary investigation has revealed that direct propagation along u of the 
subset (Ci) a to (C^) a in 1 + 3 covariant terms preserves these conditions |H. 

Recently it has been shown by Maartens |5?J and Velden ]5(| that in the subcase when p = and = u a = uj a 
('irrotational dust') all constraints are preserved along u and unique solutions exist if the initial data are analytic ( see 
also Ref. j38| ). It should be possible to extend this to the non-analytic case by enlarging the FOSH set to include the 
evolution equations for the constraints as well. Indeed, Friedrich derives and discusses a FOSH evolution system for 
the propagation of constraints of a particular FOSH reduction of the vacuum EFE based on a combination of ADM 
3 + 1 and Bianchi field equations . We anticipate that extension to at least an irrotational perfect fluid should be 
straightforward. We explicitly demonstrate this below in the particular case of cosmological models in LRS class II. 



D. Propagation of inhomogeneities 

It is instructive to take a qualitative look at how, in cosmological models (Ai, g, u) with barotropic perfect fluid 
matters source fields, changes in the distribution of the matter are communicated from one spacetime event to a 
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nearby one. For this reason let us define the 1 + 3 covariant gradient variables X a := D a u and Z a :— D a 8. It is then 
fairly straightforward to derive the following set of evolution equations (see, e.g., Refs. and jlTj ) 



X {a) = ~\QX a - a%X b 



+ e abc iu b X c -( i i+ P ) Z a 



(3.38) 



Z(a) _ c 2 D Q (D b u fc ) = (^D«[ (u b ii b ) - (a b c a\) + 2{u b uj b ) } - c 2 [ 9 Z a + a% Z b - e abc uj„ Z c + \X a ] 



{D b ii b ) + {u b u") 



■M 



3 w 



(a b c a c b ) + 2(Lu b cu b ) + n 



(3.39) 



(D a u a )Y-c 2 s D a Z a = -[ 



cPp 



(jjl +p) - c 2 + § ] 6 (D a u a ) - 2 a ab [ D< a ?i b > + u b *> } + 2 u a e abc D b u c 
-2 d2 V 



fl^+p)Qu« DaC ; 2 - 2 [ c ; 2 ^ (p + p) - cl + I ] Z" 



(/x +p) + (c; 2 - 1) ] c7 2 (M +P) e (« a «») + (c 2 + |)6 (u a u a ) ■ (3.40) 



As noticed above from Eq. ( 3.27 ) for the spacetime metric g and the matter energy density /1, we now find that also 
the (magnitude of the locally) spatial gradient of /i propagates along the integral curves of u, that is, growth or decay 
of X a is entirely determined locally. 

Where then do the sound waves occur in this first-order formulation of the equations? In fact, information on inho- 
mogeneities in the matter distribution is propagated through a dynamical coupling (determined by the local equation 
of state) of the spatial gradient of the fluid expansion to the spatial divergence of the fluid acceleration, enabling waves 
travelling at the speed of sound which then pass on information to fu rther fluid kinematical variables and, hence, by 
the kinematical interaction terms on the right-hand side of Eq. ( |3.38 ), to the density. When c s = 0, information can 
only be propagated through the Weyl curvature. Subcases for which even this vehicle of communication is unavailable 
(due to dynamical restrictions on the initial data) have been called 'silent' cosmological models f^^of ; in them, all 
initial data are merely evolved along u. 



E. Implications of the v = ± | Weyl curvature characteristic eigenflelds 

The key issue now is what to make of the characteristics along which the longitudinal components of the Weyl 
curvature travel at speed \v\ = \ relative to observers comoving with u, the latter being physically and geometrically 
uniquely defined in the perfect fluid case we consider here. They appear to represent the possibility of propagation of 
longitudinal gravitational waves at (in the given frame) one half the speed of light, thus corresponding in some sense 
to massive rather than massless particles. Is this really so, or is there some way to avoid this conclusion? We consider 
here a series of relevant issues; our conclusion is that none of them is at present sufficient to exclude such longitudinal 
gravitational waves. 



1. Can they be gauged away? 

The Weyl curvature itself is of course an invariant geometrical quantity (and is gauge-invariant for linearised 
gravitational perturbations of Minkowski spacetime or FLRW cosmologies). However, various Weyl curvature frame 
components can be set to zero if one is allowed arbitrary choice of the reference frame — see, for example, the 
discussions on the Petrov classification of the Weyl curvature of (vacuum) gravitational fields in the null tetrad 
literature, in which one of the normal forms for the algebraically general case sets the longitudinal components to zero 
fjBf , but another instead sets the transverse components to zero [0. But that frame freedom is not allowed in the 
situation considered here, where { e a } is tied in to the Fermi-propagated reference frame of fundamental observers 
- which is the relevant frame to consider in the context of detection of gravitational radiation via the GDE 
Thus, as emphasized in the previous section, one can use the remaining frame freedom to set one of the longitudinal 
polarisation modes to zero, but not to remove all components of a propagating longitudinal gravitational wave. There 
is not enough frame freedom remaining. 

Furthermore, these modes do indeed cause physically detectable effects, via the GDE (cf. the discussion in the 
introduction, and the display of the GDE source terms due to the longitudinal waves that is given in appendix |A 5| ). 
In principle, one could cancel these effects by accelerating and rotating the detector in precisely the right way, if 
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one knew in advance what information was going to be contained in the wavep| — but that does not correspond to 
motion of a physically realistic detector; and if we were somehow to achieve this, the information contained in these 
wave modes would then be coded in the requisite motion of the detector (which could only be achieved by use of 
suitable rocket motors or other energy-consuming physical actuators). Information (and energy) would indeed have 
been conveyed in a physically detectable fashion from the source to the observer. 

2. Does some consistency condition prevent them occurring? 
We have been unable to see any reason why solution of the initial constraint equations should prevent such modes 



appearing. Taking linear combinations of Eqs. (2.10) and (2.11) one finds that the propagation direction adapted, 
divergence-like ei frame derivative of the transverse modes is unconstrained whereas that of the longitudinal modes 
has typically non-zero sources. Further, we cannot see why consistency conditions between the time evolution and 
constraint equations should forbid them. To explore this a bit further, we have looked at the form of th e equations 



when linearised about a FLRW model; the Weyl curvature variables for this case are given in appendix A 4, the related 



exact propagation equations in appendix A 2. There seems to be no barrier to the linearised equations having such 
solutions. 

3. Does some conservation law prevent their emission? 

Overall conservation of mass and momentum, resulting in asymptotic conservation laws in asymptotically flat space- 
times ( see, e.g., Bondi et al || ), are global conservation laws over the entire sky. It is not clear how they could prevent 
emission of radiation of the type identified here in specific directions, provided the overall balance over the whole sky 
was achieved. 

To properly examine conservation laws associated with this kind of radiation, one would, however, have to rework 
the previous asymptotic analysis in terms of coordinates adapted to the semi-null cone corresponding to a speed of 
propagation \v\ = ^, rather than the Lorentz-invariant null cone itself. We have not attempted that extension of 
previous work, but see no reason to believe it would lead to conservation laws that would forbid radiation such as 
identified here. 

4. Why have they not been detected in previous work? 

Given that these longitudinal modes seem to be unequivocally implied by the FOSH form of the 1 + 3 covariant 
equations, why have they not been identified before? Particularly, for example, why do they not appear in the 



somewhat similar analysis of Ehlers et al 14 1, examining propagation of linearised perturbations of a 'dust' spacetime 
geometry? 

The answer seems to be that they have not been looked for, in that in most such studies the characteristics of 
propagation of the Weyl curvature were not examined; rather, those of the metric and its first derivatives were 
the focus of attention. In that case, one will not necessarily detect these modes, for the metric itself has as its 
characteristics the timelike coordinate lines associated with the fundamental observers ( cf. also Ref. ||2| ) . The key 
point is that, given a second-order hyperbolic equation for some field, the 3-surfaces across which discontinuities of 
the field and/or its first derivatives can occur are not necessarily the same as those across which discontinuities of its 
second derivatives can occur. It is the latter that we have been investigating here, because that is the level at which 
invariant physical effects occur. 

The situation may actually be analagous to that pointed out above for the matter fluid energy density, which does 
not itself immediately manifest the sound cone characteristics; these affect the density through first being manifest 
in other variables whic h th en affect the density itself ( by appearing in the kinematical interaction terms on the 
right-hand side of Eq. ( |3.3S| ) ). 



Indeed, when attention has been placed on the Weyl curvature itself, these modes have been discovered before by 
Friedrich, who examined FOSH evolution systems for vacuum gravitational fields that employed the Bianchi field 



13 Charles Hellaby (private communication) suggested a nice analogy: in principle one can shake one's head in such a way as 
to not hear a sound wave. 
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equations [|4j.[_j However, he discounted them because the related characteristic 3-surfaces of propagation were not 
Lorentz-invariant, and because an invariantly defined timelike vector field did not exist As indicated above, we do 
not believe this is a good enough reason in our barotropic perfect fluid case to claim these are not physical modes 
that can transport energy from one place to another. Our emphasis on the physics associated with the GDE indicates 
the contrary. As already remarked by Szekeres 48 , 49] , the longitudinal wave component of the Weyl curvature does 
contribute to the GDE force ter m, to give measurable effects. The relevant GDE terms are given for both null and 
timelike congruences in appendix A 5 . 

Interestingly, special propagation processes of disturbances in gravitational fields along non-null curves were de- 
scribed before by Araujo in the context of non-linear interactions between three gravitational waves in a vacuum 
spacetime geometry H. 



5. How can they have a propagation speed that is not Lorentz-invariant? 



In the specific case we consider, Lorentz invariance is broken by the existence of a cosmological fluid with a preferred 
4- velocity vector field u. Many of the variables we have identified, including the Coulomb-like components of the Weyl 
curvature, have the associated timelike curves with v = as (degenerate) characteristics; and these are certainly not 
Lorentz-invariant, but we are not surprised by their existence. Similarly, also the v = ± \ modes are not Lorentz- 
invariant. When boosting to a reference frame that does not comove with u, in the perfect fluid case the longitudinal 
Weyl curvature modes typically get mixed up with the matter fluid variables /i and p and the characteristic structure 
as well as all propagation speeds other than \v\ = 1 will look different in this framc.^| However, we determine 
propagation velocities with respect to a specific set of physical observers, and this breaks Lorentz invariance. (For 
complete ness, we give the independent Lorentz-invariant scalars one can construct from the Weyl curvature tensor in 
appendix |A gj .) 

While we are surprised by what is indicated by the equations, the non-Lorentz- invariant nature of the characteristics 
for the longitudinal modes does not by itself mean they do not occur. On the contrary, if the standard theory of FOSH 
evolution systems is indeed reliable, these are characteristics along which discontinuities and waves can propagate. 
The possibility of existence of associated physical waves has to be taken seriously unless some reason is found why 
they cannot exist. As indicated above, we are unaware of any such reason. 



IV. WORKED EXAMPLE: COSMOLOGICAL MODELS IN LRS CLASS II 



After derivation of a 44-D FOSH evolution system for generic cosmological models (M, g, u) with barotropic 
perfect fluid matter sources, based on geometrical fields defined in the 1 + 3 ONF formalism, we now turn to discuss 
in some detail a relatively simple, lower-dimensional, subcase, namely cosmological models in Locally Rotationally 
Symmetric ('LRS') class II [^6| , |l9| . This class comprises cosmologies of either spherical, plane or hyperbolic spacetime 
symmetry that are typically spatially inhomogeneous. The first to present a FOSH formulation for the spherically 
symmetrical case were Kind and Ehlers ]33]| . Their work, which, besides the matter energy density, employed dynam- 
ical field variables derived from metric functions and their derivatives, led to a 8-D FOSH evolution system subject 
to five constraints, the dependent field variables being (in their notation) u A = [ X, Y, Q, A, C, w, oj, fi ] T . The 
dynamical formulation to follow yields a result of similar structure. It represents a 8-D invariant subspace of our 
full 44-D FOSH evolution system which employs geometrical variables. Ultimately, we will explicitly include the 
conservation of the constraint equations to build a larger FOSH evolution system. This latter aspect will lead to a 
natural solution for the threading lapse function M (which we do not assume to be the constant Mo of the previous 
section). Besides these issues, our formulation shows clearly that the intrinsic geometry of the 2-D spacelike symmetry 
group orbits can be fully decoupled from the evolution of the physically relevant field variables. Due to the spacetime 
symmetry imposed, these models exclude all (forms of) gravitational radiation; thus, apart from sound waves, these 
are gravitationally 'silent' universes, as mentioned above. 



1 He does not mention these modes in his fluid paper [|25[; but they must be there in the system he has examined in that 
paper, for they were there in the vacuum case. 
15 Helmut Friedrich (private communication). 

16 The given v = and v = ±c s modes are not Lorentz-invariant either; however, in the generic case there are, e.g., Coulomb- 
like components of the Weyl curvature with \v\ = in any rest frame. 
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In the 1 + 3 ONF formalism, LRS spacetime geometries are characterised as follows pi] : the unit tangent of 
the invariantly defined preferred spacelike direction may be taken as, e.g., e = ei; then = e2(/) = e 3 (/) for 
any invariantly defined geometrical variable. Rank 2 symmetric tracefree te nsors A orthogonal to u specialise to 
= A- = A x = A 2 = A 3 => A 2 = | (A + ) 2 => A = ±^A + (cf. Eq. (^2?|)). The case of our interest, the 

cosmological models in LRS class II, are defined by the conditions = u> = nu => H + = [[l6 21 1. Then it can be 
shown that the general relations listed in section [n] reduce to ( cf. Ref. [jl7| ) 



The commutators: 



The evolution equations: 



[e ,ei] = iie - ±(9 - 2<r + )e 1 

[e ,e 2 ] = - |(9 + (7+) e 2 

[e ,e 3 ] =-i(e + fr + )e 3 

[ei, e 2 ] = ae 2 

[e 2 ,e 3 ] = 2n 3 ie 3 

[e 3 , ei ] = - ae 3 . 



e (9) = - | 9 2 + (ex + u - 2a) (u) - § (a+ f - § (fj, + 3p) 
e (o+) = - | (26 - cr + ) (T + - (ex + ii + a) (u) - E + 



e (u) = (c 2 s - |) Qu + c 2 s ex (6) + © || eifc) + - 

e (a) = - ^(6 + (7+) {u + a) 
eo(n3i) = - |(S + cr + ) n 31 
e (£+) = - § (/x + p) o-+ - (9 + (7+) E+ 

eo(^) = -eo + p) . 



I (7+ ?i 



The constraint equations: 



= 


(Cr) := 


(ea 


- 3a) (*+) ^ 


-ex(9) 


= 


(C 2 ) := 


(ei 


- a) (a) + | 


9 2 - i (9 + 2o- + ) o-+ + i£; 


= 


(Cs) := 


(ei 


~ a) (n 31 ) 




= 


(C 4 ) := 


(ei 


- 3a) (£ + ) - 




= 


(Cpf) : 


= c 2 s 


ei(M) + (M- 


fp) « . 



(4.1) 
(4.2) 
(4.3) 
(4.4) 
(4.5) 
(4.6) 



(4.7) 
(4.8) 

(4.9) 

(4.10) 
(4.11) 
(4.12) 
(4.13) 



(4.14) 
(4.15) 
(4.16) 
(4.17) 
(4.18) 



Note that Eqs. ( 4.11 ) and ( 4.16| ) for the spatial commutation function n 31 decouple from the remaining set of 
equations; n 3 i in fact determines the 



Gauflian curvature of 2-D spacelike symmetry group orbits: 

K =2 (ea - 2n 3 i) (n 3 i) ; 



(4.19) 



the group orbits can either be spherically symmetrical (K > 0), plane symmetrical {K = 0), or hyperbolically 
symmetrical (K < 0). 



Tracefree part and trace of 3-Ricci curvature: 

*S+ = - ei(o) + K = E+- |(9 + a+) a+ + (C G )+ 
*R=2(2e 1 -3a)(a) + 2K = 2 fi - § 9 2 + § (<r+) 2 + (C G ) 



Combining Eqs. (4.15) and (4.20), one obtains for K the algebraic expression 

17 ' " |(9 + or+f + a 2 + (C G )+ + (C 2 ) 



(4.20) 
(4.21) 



(4.22) 
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A. Set of FOSH evolution equations 



We choose our set of (dimensionless) comoving lo cal coordinates as { x^ } = { t, x, y, z }. Since cu a = => Mi = 0, 
the expansion of the frame basis fields e a of Eq. (3.6) reduces to a zero-shift-vector version ('Eulerian observers') 
of the ADM 3 + 1 slicing formulation of the dynamics of relativistic spacetime geometries [Q ; with properly defined 
spacelike 3-surfaces T: {t — const}. For LRS class II we find in particular 0,^6) 



e = M" 1 d t , e 2 = Y- 1 d y , e x = A" 1 d x , e 3 = (Y Zy 1 8 Z , 



(4.23) 



where M — M(t, x), X = X(t, x), Y — Y(t, x) and Z = Z(y); each of the e a here is hypersurface orthogonal. With the 
simplifications due to the spacetime symmetry of LRS class II given, our set of dynamical field variables is comprised 
of the 8-D subset u A = [ii, 1(8 - 2cr + ), |(8 + o+), a, E+, fi, A, Y } T of Eq. ( j3.17fl , leading to the FOSH evolution 
system 



M" 1 8 t u - c 2 A" 1 94(9 - 2a+) = - [ c; 2 ^| (/x +p) - c 2 + 



V ■ *v -» 1 3 



9; 



c 2 (d 



8 c: 



d 2 P 
4? 



- ^a + u + 2c s acr + 
(C PF ) 



c 2 M" 1 9 t i(6 - 2cr+) - c 2 A" 1 c^u : 



c 2 [ i (6 - 2a + ) 2 - u 2 + I (/i + 3p) - § £ + 



UQ + a+){u + a) 



M- 1 5 t i(9 + o+) = - i (6 + a+f - aii-\{ji + Zp)-\E + 
M~ l d t a ■■ 
M~ 1 d t E + : 



±(m + p) ( T + -(e + ( 7 + )s+ 
e(/i + p) 



(4.24) 
(4.25) 
(4.26) 
(4.27) 
(4.28) 
(4.29) 



M^d t X = 1(6 - 2<j + )X 



(4.30) 



M -1 sty = i(e + cr + )r 



(4.31) 



We note in passing that, like in Eq. ( 4.11 ), the evolution of the variable Y, which will be found to be equivalent 
to the Gaufiian curvature according to K — C\/Y 2 , is decoupled from the remaining equations; thus the resulting set 
is effectively 7-D. 

The symmetric matrices M AB ^ entering in Eq. (1.6) assume in the given situation the explicit forms 



/I 



M 



ABO 



( ■ 1 
1 . 



M AB1 = -MX~ 1 ct 



1/ 



V 



(4.32) 



/ 



The FOSH evolution system ( 4.24 ) - ( 1.31 ) will be fully deterministic, if either wc prescribe the threading lapse 
function M by hand, or we have a further equation for its evolution along u. We will return to this issue further 
below. 

The characteristic condition = det [ M AB M £ M ] , with, in the present case, ^=dM 5°^ + X S 1 ^, reads 



The characteristic condition: 

= det [ v M M AB + X M AB 1 ] = M 8 c 2 v 6 (v - c.) (v + c s ) 



(4.33) 
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Hence, the characteristic propagation velocities with respect to a n obser ver at rest in { e Q } are either v\ — or 
U2.3 = ± c s . Upon diagonalisation of the principle part of Eqs. ( 4.24 ) - ( 4.31 ), the characteristic eigenfields propagating 
along the sound cone are given by 



'(2) 



y/T+c 



(c s u±U®-2a + )) ] 3 



(4.34) 



a res ult w hich, as to be expected, is just a specialisation of the general one obtained for barotropic perfect fluids in 
Eq. (|3.28|) of section |ffl]. 

Setting p = 0=>« = 0, c s =0 (for K > 0), we get from Eqs. (4.24) - (4.31) the evolution part that describes the 
('silent', and spherically symmetrical) Lemaitre Tolman Bondi cosmological models. 



B. Propagation of constraints 

Using the coordinate components of { e a } given in Eq. ( [4.23| ), the commutator relations ( f4.1| ), ( |4.4| ) and (4.5) yield 



0=(Cy := u- {MX)~ 1 d x N 
0=(C„) := a+iYX^dxY 
0=(C„ 31 ) := n 3 i + \ (ZY)- 1 d y Z , 



(4.35) 
(4.36) 
(4.37) 



respectively. Together with Eq. ( 4.37| ), Eq. ( 4.1 9| ) integrates to K = Ci/Y 2 , where Ci deno tes a dime nsi onles s 
integration constant. One can show that the set of constraints comprised of Eqs. (4.14) - (4.18) and (4.35) - (4.37) 
propagates according to 



M- 1 d t (d) = -Q (Ci) - 3 u (C 2 ) - (C 4 ) - § (C PF ) 



+ [ i 8 2 + § Sa+ + \{a+f + 3au+ \ (jm + 3p) + E+ } (C fi ) 



i (29 - a+) (C 2 ) -l(u + a) (d) + 1(6 + a + ) (u + a) (C*) 
i (29 - a+) (C 3 ) - | n 31 (Ci) + |(6 + n 31 (C 4 ) 



Af- 1 3 t (C 4 ) = - i (46 + o-+) (C 4 ) - | (/x + p) (Ci) - £+ (Ci) - | (9 + a+) (C PF ) 

+ \ (9 + <r+) (a* + V) (Cu) + (6 + 0-+) E+ (Cu) 

M- 1 d t (C PF ) = - [ (c 2 + |) 6 + § a + } (C PF ) + c 2 s 6 (/x + p) (C 4 ) 

M" 1 9 t (d) = - |(6 - 2a+) (C,0 - [ (MX)" 1 d t {M- x d x M) - X- 1 d x (c 2 s e) ~c 2 s Qu] 

M- 1 d t (C a ) = - |(6 - 2a+) (C a ) + | (Ci) - |(6 + a+) {Cu) 
M- 1 d t (C n31 ) = -Ue + a + )(C n31 ) . 



(4.38) 
(4.39) 
(4.40) 

(4.41) 
(4.42) 
(4.43) 
(4.44) 
(4.45) 



The constraints will propagate consistently, in fact also via a FOSH evolution system (with zero characteristic prop- 
agation velocites only), iff the evolution of the threading lapse function M (i.e., d t M) satisfies the relation 



OtlM^dxM] = Md x [c 2 e] + MXc 2 Qu , 



(4.46) 



which follows from Eq. (4.43). That is, the evolution of M is de term ined through an integrability conditi on de manded 
by preservation of the constraint s. It t urns out that th e cho ice ( |3.15 ) of Salzman and Taub satisfies Eq. ( 4. 46] ) subject 
to the validity of the constraint ( 4.35| ). The solution (3.15) implies 



NT 1 d t M = c 2 QM 



(4.47) 



Given the results of this subsection, it is natural to add Eqs. ( 4.38| ) - (4.45) to the set (4.24) - ( [4.31 ) to form a 
larger FOSH evolution system for the set of dynamical fields 
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^ A = t yfe ( Cs " ± ^ e - 2<7 +) )i 3( + CT +)' a > S +> M. X, Y, (C n ) ] T , (4.48) 
which can be taken to model all spacetime geometries with perfect fluid matter source fields in LRS class II. 

V. CONCLUSION 

The aim of the present work is to explore the propagational features of physically relevant quantities in relativistic 
cosmological models with idealised matter source fields that can be treated as perfect fluids with barotropic equations 
of state. The results we obtain on the characteristic structure of the derived 44-D FOSH evolution system apply to 
both, exact cases and FLRW- linearised cosmological models; the principle part (the dynamical interactions between 
various fields) remains unchanged under the linearisation procedure, which affects the kinematical interactions only. 
Hence, we can read off which kind of perturbations propagate along which set of characteristic 3-surfaces in both the 
exact and the linearised theories. 

As to the question raised in Ref. |2Q] concerning the presence of gravitational radiation in spatially homogeneous but 
anisotropic cosmological models with perfect fluid matter source fields, it now becomes clear that, with respect to a 
symmetry-group-invariant 1 + 3 ONF { e a }, there then do not exist any associated null characteristics (while they are 
allowed by the dynamical equations, the homogeneous data does not activate them). Only timelike ones are present, 
each with propagation velocity v = of the associated characteristic eigenfields. Hence, there exists no propagating 
gravitational radiation in such models (nor any other fields that are exchanged between neighbouring worldlines of 
the matter fluid). In this sence, this class of cosmological models can be classified as 'silent', too. For certain cases 
in this class it may, however, be sensible to maintain some kind of interpretation of the spacetime dynamics in terms 
of standing gravitational waves (with respect to the timclike reference congruence chosen); aspects along these lines 
have been discussed by, e.g., Lukash or King Q. 



The 1 + 3 ONF formulation with its geometrically defined dependent field variables provides a possible alternative 
FOSH evolution system for applications in numerical relativity, in particular to relativistic cosmology ( for a 'living 
review' on the current state of the latter see Ref. ). Even though, in this framework, the number of constraints, 
which are typically non-linear elliptic partial differential equations, is larger than in conventional schemes based on 
the ADM 3 + 1 treatment, the benefits may be substantial. Most variables in our treatment have direct physical 
significance (e.g., the Weyl curvature and matter variables enter the GDE); hence, they are intuitively more easily 
accessible than the derivatives of metric components which appear in various other works. The example of cosmological 



models in LRS class II described in section IV may serve as a simple test bed for numerical simulations; analysis of 
'shell-crossing singularities' in the Lemaitre-Tolman-Bondi 'dust' case being one of possible problems to investigate. 
Moreover, irrespective of the fact that a timelike reference congruence u is not invariantly defined for vacuum spacetime 
geometries, there is no reason that prevents application of the 1 + 3 ONF FOSH evolution system to those situations 
as well, as recently emphasised by Friedrich |24| . For example, in Ref. (4f| a 16-D FOSH evolution system is presented 
for the numerical description of a static Schwarzschild black hole in harmonic time slicing. Modifying our general 



equations for /i = along the lines outlined for spherically symmetrical cases in section IV could lead to a lower 
dimensional evolution system which may, furthermore, not suffer from the continuum instabilities reported in that 
article. There the authors attribute this problem to the intrinsic structure of the equations they employed. That the 
FOSH evolution system of the 1 + 3 ONF formalism does not have flux-conserving form may or may not have negative 
consequences in prospective numerical investigations (cf. Ref. | p4| and the 'living review' p3| ). 

The most unexpected feature arising from our FOSH analysis of invariantly defined quantities is the indication of 
the existence of propagating longitudinal modes in the Weyl curvature. If this is confirmed by further analysis, it has 
possible intriguing implications in regard to (a) gravitational wave emission, (b) gravitational wave detection, and (c) 
in relation to the quantisation of the gravitational field. As regards the first two issues, they presumably provide a 
further channel whereby energy can be emitted from a source in the form of gravitational radiation and then received 
in distant places (hence, resulting in detectable effects). This would mean they need, for example, to be considered 
seriously in the context of the inflationary universe scenarios for gravitational wave emission and the associated CMBR 
anisotropics. As regards the latter, they imply a substantial difference relative to the vacuum electrodynamical case, 
where longitudinal modes exist but have no physical effect and can be gauged away. However, in a plasma, where 
(charged) matter is present as well as the electromagnetic field, the situation seems rather similar to what is derived 
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here. 17 We have shown that in the non- vacuum gravitational case longitudinal modes appear that cannot be gauged 
away; hence, possibly, these modes may need to be taken into account when quantising gravitational radiation in 
some contexts. 
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APPENDIX: 
1. Conventions 

Throughout our work we employ geometrised units which are defined by setting c = 1 = 8irG/c 2 . All geometrically 
defined variables thus have physical dimensions that are integer powers of [ length ] only. We denote spacetime indices 
with respect to an arbitrary basis { e a } by a, 6, c, ... = — 3. When fixing the basis vectors to be orthonormal 
such that the components of the metric tensor g reduce to the Minkowskian form g ab = rj ab = diag ( — 1, 1, 1, 1 ), and 
aligning eo with the future-directed, normalised tangents u of a congruence of preferred timelike curves, indices with 
respect to the spatial basis { e a } are denoted by a, j3, 7, . . . = 1 — 3. For a set of dimensionless local coordinates 
we write { x 1 * } = { t, x l }; spacetime and spatial coordinate indices being fi, v, p, . . . = — 3 and i, j, k, . . . = 1 — 3, 
respectively. 

We define the metric compatible covariant derivative V and so the spacetime connection T a b c by 

V a V b := e a (V b ) + T b ca V c , V a V b := e a (Vfc) - r c fca V c . (Al) 

In an orthonormal frame { e a } we have 

= V a .g 6c = V a % c = - T d ba j] dc - T d ca rj M , (A2) 

i.e., F( afc ) c = 0. For easier direct comparison it should be noted that Friedrich's notation for the connection components 
is V a V b := e a (V b ) + T a b c V c instead 24 25|| , such that the index referring to the 'left-acting' derivative is the first one 



to appear on r a fc c . Compared to Ref. ]21|7we will choose the fluid vorticity vector u a to have the opposite sign; thus, 
its definition will now comply with the Newtonian convention. Finally we should remark that the Weyl curvature 



variable B a b used by Friedrich [|24 25 has the opposite sign to our H a b, as in his definition the negative of the 3- volume 
element was used. Our convention ensures the same relative sign structure between E ab and H a b in the Bianchi field 
equations when compared to that between E a and H a in the Maxwell field equations. 

2. Evolution of Weyl curvature characteristic eigenfields 

Transverse v — ±1 modes: 

(e T ei) (£?_ =F Hi) ± \ e 2 {E 3 T H 2 ) T \ e 3 {E 2 ± H 3 ) = - (9 - 0+ T 2 u x ± a x ) (£_ T #1) 

- {uj 1 + 2n 1 ±n 22 ±n 33 T^n 11 ){E 1 ±H-) 

- \ (V3cr 2 - lj 2 - 2fl 2 t 2« 3 ± a 3 ±3n u ) (E 2 ± H 3 ) 
+ \ (n/3(7 3 + uj 3 + 2fl 3 T 2ii 2 ± a 2 =F 3 n 3 i) (E 3 =p H 2 ) 

- \ (n+p)a- + (cr_ =F Vin 23 )E + 



7 Malcolm MacCallum (private communication). 
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T (a 1 ±^n 22T ^n 33 )H + (A3) 

(e T ei) (E 1 ±H_)±\ e 2 (E 2 ± ff 3 ) ± I e 3 (£ 3 T^ 2 ) = -(9- ( 7 + T2 I i 1 ± oj) ± fl_) 

+ (wi + 2 n x ± n 22 ± n 33 =F | nn) T Hi) 

+ \ (VSo-i -LU 2 -2n 2 T2u 3 ±a 3 ±3n u ) (£3 =p #2) 

+ i (V3a 3 + co 3 + 2fi 3 T 2ii 2 ± a 2 T 3n 3 i) (#2 ± As) 

- +p) <7i + (cri ± ^ n 22 T n 33 ) E+ 

± (a_T\/3n 23 )ff+ . (A4) 

Note that, with the exception of the source terms proportional to er_, ci, £'+ and -H+, the outgoing characteristic 
eigeniields are driven only by outgoing characteristic eigenfields and, similarly, the incoming characteristic eigenfields 
are driven only by incoming characteristic eigenfields. 



Longitudinal v = ± ^ modes: 

(e T | ei) (£ 3 T H 2 ) ± | e 2 (£_ T H x + V3 E+) ± I e 3 (^ ± #_ ± x/3 ff+) 
= -i(20 + a + T2i 11 ± fll )(B 3 TiJ 2 ) 

+ ((7- T V3 71 23 ) (£ 3 ± fl2) - 3 (Wi + 2 fii) (£ 2 ± ff 3 ) + ^ (Ji (£2 T # 3 ) 

+ i (V3 a 3 - uo 3 - 2 n 3 T 2 u 2 ± a 2 ± 3 n 31 ) (J5L T #i) 
+ I (V3 era + cj 2 + 2 Q 2 T 2 u 3 ± a 3 T 3 m a ) (£1 ± 

- i (/j + p) ct 3 ~ (^5 ^3 + ^3 + 2 fi 3 =F 2 ii 2 ± a 2 =F «3i) ^+ 
± ^(^ f fT 2 -c i ; 2 -2r! 2 T2ij3±a3±^i2)H + 

T I (2 n 33 + 2 nn - n 22 ) £2 - \ (2 nn + 2 n 22 - n 33 ) H 3 (A5) 

(eo T § ei) (£2 ± H 3 ) ± i e 2 (£i ± # - T V3 fl+) T § e 3 (£_ T H x + V3E+) 
= - \ (2 6 + a+ t 2 uj ± 01) {E 2 ± JT 3 ) 

- ^ (er_ T V3 n 23 ) (£2 T As) + I (wi + 2 Hi) (£ 3 T i/ 2 ) + | a x (E 3 ± ff 2 ) 

- i(\/3 < 7 2 +^ 2 + 2f7 2 T2 7j 3 ±a 3 T3ni 2 )(£;_ T^i) 
+ i(V3(J3-W3 -20 3 T2« 2 ia 2 i3 n 3 i) (#i ± H-) 

- \ {ji + p) a 2 - ^ a 2 - uj 2 - 2 Q 2 T 2 ii 3 ± a 3 ± m 2 ) 

T ^ (75 CT 3 + ^3 + 2 ft 3 T 2 ii 2 ± a 2 T » 3 i) 

± I (2 nn + 2 n 22 - n 33 ) £3 - § (2 n 33 + 2n n - n 22 ) H 2 . (A6) 
3. Lorentz-invariant Weyl curvature scalars 



In terms of the tracefree-adapted frame variables of Eq. (2.27) and the characteristic eigenfields of the Weyl 



curvature, the four standard Lorentz-invariant scalars assume the forms 

1 r~iab r~icd rpa rpb rja rjb 

s° cd^ ab — ^b&a — tlbtla 



§ [ (E+f + (H+f + (E_ - H x ) (E- + Hi) + (E 1 + HJ) {E x - H-) 



+ (E 3 ~ H 2 ) (E 3 + H 2 ) + (E2 + H 3 ) (E2 - H 3 ) ] (A7) 



*t> cd<^ ab — — 1£j btl a 

_ 2 
3 



2 [ 2 E+ H+ - (E- - Hi) (Ex - H-) + (E- + H t ) (E 1 + H-) 
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- (E 3 - H 2 ) (E 2 - H 3 ) + (E 3 + H 2 ) (E 2 + H 3 ) } 



(A8) 



1 r^ab /^icd r~ief j^a jpb tjc , o pa rjb rjc 

YgO c dLy e f O ab — — ^b^c^a+OtL},tl c tl a 

= | [ i (E+f - E+ (H+) 2 - E+ (£_ - H{) (E_ + ff x ) - + - Jf_) 

+ i £+ (£ 3 - JT 2 ) (E 3 + H 2 ) + \ E+ (E 2 + H 3 ) (E 2 - H 3 ) 

- H+ (£?_ - ffx) (Bi - £T_) + (E_ + Hi) {E 1 + HJ) 
+ \H+ (E 3 - H 2 ) (E 2 -H 3 )-\ H+ (E 3 + H 2 ) (E 2 + H 3 ) 

- ^ (E_ - Hx) (E 3 + H 2 f -^(E-+Hi.) {E 3 - H 2 f 
+ &(E_- Hx) (E 2 - H 3 f + ^(E_+ Hx) {E 2 + H 3 f 

- & (Ex + HJ) (E 3 + H 2 ) (E 2 - H 3 ) 

- y£(Ex-H-)(E 3 -H 2 )(E 2 + H 3 )} (A9) 

i *C ab cd *C cd ef *C ef ab = H% H\ H c a - 3 H\ E b c E c a 

= - § [ | [H+f + H+ (E+) 2 + H+ (E_ - Hx) GEL + H x ) + H+ (Ex + HJ) (Ex - HJ) 
-\H+ (E 3 - H 2 ) (E 3 + H 2 )-\ H+ (E 2 + H 3 ) (E 2 - H 3 ) 

- E+ (E_ - Hx) (Ex -HJ) + E+ (£_ + Hi) (Ex + HJ) 
+ \E + (E 3 - H 2 ) (E 2 - H 3 ) -\E + (E 3 + H 2 ) (E 2 + H 3 ) 

+ Jf(Ex+HJ) (E 3 + H 2 f -JZ(Ex-HJ) (E 3 - H 2 f 

- ^ (Ex + HJ) (E 2 - H 3 f + ^(Ex- HJ) (E 2 + H 3 ) 2 

- & (E_ - Hx) (E 3 + H 2 ) (E 2 - H 3 ) 

+ & (E_ + Hx) (E 3 - H 2 ) (E 2 + H 3 ) ] . (A10) 



4. Weyl curvature variables in linear perturbation formalisms of FLRW 

Linearised perturbations to FLRW cosmological models are often analysed in terms of the non-geometrical gauge- 
invariant variables suggested by Bardeen || . From the paper by Goode (3^] ( see also Ref. [|l0| ) it is straightforward 
to expand the characteristic eigenfields of the Weyl curvature in terms of the former. In the absence of dissipative 
perturbations to the perfect fluid matter content it follows that $,4 = — . 

Scalar modes: 



(All) 



Vector modes: 



E+ 


= -(3/4)fc 2 ($ A ™ 




£L 


= (V3/4) fc 2 ($ A - 


$ ff )[Q (0) 22-Q (0) 33 


Ex 


= (V3/4) fc 2 ($ A - 


$H)Q {0) 2 3 


E 2 


= (V3/4) fc 2 ($ A - 


$H)Q (0) 31 


E 3 


= (V3/4) fc 2 ($ A - 


$H)Q (0) 12 


H a 


= . 





E+ = -(3/4)fc#Q« u 

H+ = -(3/4)*[QW 2 | 13 -QW3| 12 ] 

(E.tHx) = (V3/4)fc*[Q(i) 22 _Q(i) 33 ] T (V3/4)*[g( 1 ) 1 | 22 -Q( 1 ) 1 |33-Q (1) 2 |2i+Q (1) 3|3i ] rAl2 ^ 

(Ex±HJ) = (V3/2)fcvi/Q( 1 ) 2 3±(v / 3/4)*[g( 1 ) 2 | 31 + Q( 1 )3|2i-Q (1) i|23-Q (1) i|32] 1 1 

(E 3 T H 2 ) = (V3/2) k * QW 12 T (73/4) * [ QW 2 | 33 - Q( 1 ) 2 | 11 - Q( 1 ) 3 | 32 + Q (1) i|i2 ] 

(E 2 ± H 3 ) = (V3/2) k * Q« 31 ± (^3/4) * [ Ql» m - Q^ 3 \ 22 - Q^ l \\ 13 + Q^ 2 \ 23 ] . 
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(A13) 



Tensor modes: 

E+ = (3/4) [ - (k 2 + 2K) H {2) } u 

H+ = - (3/2) [ Q^) 12 | 3 - Q (2) 3i|2 ] 

(E_ T Hi) = - (V3/4) [ H ( t ] - (k 2 + 2K) H {2) ] [ Q^ 22 - Q^ 33 } 

± (V3/2) [ Q^ 22]1 - Q( 2 ) 33 , 1 + Q( 2 ) 31 | 3 - Q( 2 ) 12 , 2 ] 
(E 1 ±H_) = - (V3/4) [ - (k 2 + 2K) i4 2) ] g( 2 ) 23 

± (V3/2)^ 2) [ 2Q( 2 ) 23 , 1 - Q( 2 ) 31 , 2 - Q( 2 ) 12 , 3 ] 
(E 3 =F H 2 ) = - (V3/4) [ HP - (k 2 + 2K) ] Q^ 12 

± (V3/2) [ Q( 2 ) 33 , 2 - Q( 2 ) U | 2 + Q^ 12]1 - Q( 2 ) 23 , 3 ] 
(E 2 ±H 3 ) = -(V3/4)[^ 2) -(fc 2 + 2X)4 2) ]g( 2 ) 31 

T (V3/2) ^ [ Q^llla - Q (2) 22|3 + (2) 23|2 - (2) 31|1 ] • 

5. The Weyl curvature force term in the deviation equation for null and timelike congruences 

We choose the tangents to the null curves to be given by k a := E (u a + e a ), where E := —(u a T] a ). Then the 
orthogonal deviation vector can be expressed as rf = — E^ 1 (ubi] b ) k a + r\ 2 e 2 a + n 3 e 3 a . Thus, the force term is 

Null congruences: 

C a bcd k b n c k d = ^E{ (E 3 - H 2 ) m + (E 2 + H 3 ) n 3 } k a + ^ E 2 [ (E_ - Hj) m + (E l +H_) m ] e 2 a 

+ ^=E 2 {(E 1 + H_) m - (E_ - H{) m } e 3 a . (A14) 

For timelike curves moving with velocity v a :— ve a relative to u we define the unit tangents by V a := 7 (u a + ve a ), 
where 7 := (1 — tj 2 ) -1 / 2 . Then the orthogonal deviation vector can be expressed as rf = {e b n b ) (vu a + e a ) + n 2 e 2 a + 
i] 3 e 3 a . Thus, the force term is 

Timelike congruences: 

C a bcd V b V c V d = -j={j= E+ r ll - 1 2 (E 3 -vH 2 ) m - 7 2 (E 2 + v H 3 ) m }(vu a + e a ) 

+ [ (E 3 - v H 2 ) Vl + 7 2 (-^ E+ - v 2 E+ + E_ + v 2 E_ - 2 v H x ) n 2 

+ 7 2 {E 1 +v 2 E 1 + 2vH_)r l3 } e 2 a 
+ ^[(E 2 + vH 3 ) m +-f 2 (E 1 +v 2 E 1 +2vH-)r l2 

+ 7 2 E+ - ^ v 2 E+ - E_ - v 2 E_ + 2 v H{) r, 3 } e 3 a . (A15) 
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